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Abstract 

Exactly solvable models of planar polygons, weighted by perimeter and area, have deep- 
ened our understanding of the critical behaviour of polygon models in recent years. Based 
on these results, we derive a conjecture for the exact form of the critical scaling function for 
planar self-avoiding polygons. The validity of this conjecture was recently tested numerically 
using exact enumeration data for small values of the perimeter on the square and triangular 
lattices. We have substantially extended these enumerations and also enumerated polygons 
on the hexagonal lattice. We also performed Monte-Carlo simulations of the model on the 
square lattice. Our analysis supports the conjecture that the scaling function is given by the 
logarithm of an Airy function. 

1 Introduction 

The model of self-avoiding polygons (SAPs) is an important unsolved combinatorial problem in 
statistical mechanics [241 I13| . In this article, we consider the model of SAPs in two dimensions. 
SAPs counted by perimeter are the canonical model of ring polymers, while when also counted 
by area, they serve as a model of vesicles The model also has connections with many other 
models of statistical mechanics where boundaries of certain cluster types play a role, notably in 
percolation and in spin models such as the Ising model and the Potts model, see also [7j. Despite 
its importance, there are only a few known rigorous results about SAPs, almost all information 
arising from numerical investigations ^2 1241 1111 1131 |H1 E] • Promising new developments however 
arise from the theory of stochastic processes U^J, see also |S] (this volume). Recently, we verified 
numerically a conjecture of one of us about the scaling function for SAPs, which describes the 
singular part of the SAP perimeter and area generating function about its tricritical point |36U37| . 
The prediction of the scaling function also gives the distribution of area moments. The conjecture 
was already stated in 1995 |32| . as a result of a systematic investigation of simple exactly solvable 
subclasses of SAPs. There it was noticed that the critical exponents of rooted SAPs coincide with 
the exponents of staircase polygons, where the scaling function had been computed |29|. This led 
to the question whether both scaling functions are of the same type. 

Over the years, it turned out that the mathematical structure underlying exact solvability 
is given by g-algebraic functional equations. They arise from the requirement that the polygon 
model is built up recursively, see e.g. [Q. This general type of functional equation appears 
first, within the framework of algebraic languages, in [5], together with the discussion of the area 
distribution for simple polygon models showing a square-root singularity. However, the article |5] 
does not contain a derivation of its central result. Recently, we investigated the scaling behaviour 
of general g-algebraic functional equations [3^1 E|- The techniques used for handling g-algebraic 
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functional equations are applicable in a wide variety of problems, see e.g. |18j for exact low 
order cluster moments of asymmetric compact directed percolation. An important outcome of our 
previous analysis was that the scaling behaviour generically depends only on the singularity of the 
perimeter generating function. This remarkable kind of universality further supported the above 
conjecture and led us to test it numerically. 

Very recently, investigations on the self-avoiding walk problem using the theory of stochastic 
processes led to a number of predictions of the behaviour of SAPs in the scaling limit. One of them 
is that the scaling limit of the half plane infinite SAP is the outer boundary of the union of two 
Brownian excursions from to oo in the upper half plane |19|. Since the area under a Brownian 
excursion is Airy distributed |2U, this would suggest that the area moments of SAPs obey an 
Airy type distribution. This is consistent with our conjecture that the SAP scaling function is 
the logarithm of an Airy function and also illuminates the connection with the model of staircase 
polygons, whose scaling limit is described by Brownian excursions. 

In the remainder of this article we first introduce polygon models, in particular self-avoiding 
polygons, and discuss their critical behaviour. We then consider g-algebraic functional equations 
and show how to extract scaling behaviour in the vicinity of a square-root singularity from the 
functional equation. The following section analyses the exactly solvable example of bar-graph 
polygons, using the techniques described in the previous section. The next section on self-avoiding 
polygons contains a summary of the analysis of recently obtained exact enumeration data for 
polygons on the square lattice to perimeter 100, the hexagonal lattice to perimeter 140 and the 
triangular lattice to perimeter 58. The analysis confirms the conjectured form of the scaling 
function in that predicted exact amplitude combinations are confirmed to within a numerical 
accuracy of 5-6 significant digits. This section also explains the Monte-Carlo simulation method, 
which we performed on the square lattice for polygons up to perimeter 2048, and confirms the 
conjectured form of the scaling function within a numerical accuracy of 3 significant digits. A 
concluding section discusses possible applications of our results and methods. 

2 Polygon models and their scaling behaviour 

For concreteness, we discuss models of polygons defined on the square lattice 1? . These models 
can be defined similarly for other lattices as well as in the continuous case. A self-avoiding polygon 
(SAP) is a closed, non-intersecting loop on the edges of the square lattice. Important subclasses 
are staircase polygons and bar-graph polygons j.'iOj , whose obvious definition may be extracted from 
Figure H Whereas the latter two models have been solved exactly, as also discussed below, for 
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Figure 1: Prominent polygon models. From left to right: self-avoiding polygons, staircase polygons, 
bar-graph polygons. 

SAPs there are very few rigorous results, complemented by a number of numerical investigations. 
Let p m ,n denote the number of different polygons of perimeter m and area n, respectively, where 
we regard translated polygons as equal. The object of interest is the perimeter and area generating 
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function 

G(x,q)=Y / Pm, n x m q n , (2.1) 

where we introduced a perimeter activity x and an area activity q. Sometimes, horizontal and 
vertical perimeter are distinguished, leading to an anisotropic perimeter and area generating func- 
tion G(x, y, q). The function G(x, 1) is called perimeter generating function. Of particular interest 
is the singular behaviour of G(x, q) as a function of x and q, which translates to the behaviour of 
the coefficients p m .n in the limit of large perimeter and area. 

The above models display a so-called collapse transition between a phase of extended and a 
phase of deflated polygons 11 : in the extended phase q = 1, the mean area of polygons (a) m of 
perimeter m grows asymptotically like m 3 / 2 , whereas it grows like m in the deflated phase q < 1. 
It can be shown that in the limit q — > the generating function is dominated by polygons of 
minimal area. Since for SAPs these polygons may be viewed as branched polymers, the phase 
q < 1 is also referred to as the branched polymer phase. This change of asymptotic behaviour 
is reflected in the singular behaviour of the perimeter and area generating function. Typically, 
the line q — 1 is a line of finite essential singularities for x < x c . The line x c (q), where G(x, g) 
is singular for q < 1, is typically a line of algebraic (staircase, bar-graph) or logarithmic (SAPs) 
singularity. For branched polymers in the continuum limit, the logarithmic singularity has been 
recently proved in j^] (see also this volume). 

Of special interest is the point (x c , 1) where these two lines of singularities meet. The behaviour 
of the singular part of the perimeter and area generating function about (x c , 1) is expected to take 
the special form 

G(x, q) ~ G<™\x, q) + (1 - q) 9 F((x c - x)(l - q)~% (x, q) -> (x~ , 1"), (2.2) 

where F(s) is a scaling function of combined argument s — (x c — x)(l — q)~^, commonly assumed 
to be regular at the origin, and 9 and 4> are critical exponents. The singular behaviour about q = 1 
at the critical point x c is then given by G^ sln9 \x c , 1) ~ (1 — q) e F(0). This scaling assumption 
implies an asymptotic expansion of the scaling function of the form 

OO j. 

^E^fSw (2-3) 

k=Q 

The leading asymptotic behaviour characterises the singularity of the perimeter generating func- 
tion via G(x, 1) ~ fo(x c — x) -7 , where + 07 = 0. The first singularity of F(s) on the negative axis 
determines the singularity along the curve x c (q). The locus on the axis (say at s = s c ) determines 
the line x c (q) ~ x c — s c (l — q)^ near q — 1, which meets the line q — 1 vertically for <f> < 1. 

The area moments gk (x) are defined as coefficients of the (formal) expansion of the generating 
function about q = 1 

G(x, (? )=^ 5fe (x)(l-g) fc , g k (x) = ^^n(n-l) fc + l)p m ,„x m . (2.4) 

k ' m,n 

Note that go(x) = G(x, 1) is the perimeter generating function. The higher moments <7fc(x) can 
be computed recursively using gi(x), where I < k using the functional equation There is 

an important relation between the scaling function and the area moments: The coefficients fa in 
the asymptotic expansion of the scaling function are the leading singular amplitudes of the area 
moments g k {x). This follows from the scaling assumption which implies that the leading singular 
behaviour of the area moments is of the form 

9 { r 9 \*)~J^^ k (*-x-), (2.5) 

where jf. = (k — 9)/<j>. Thus, information about the scaling function can be obtained from the 
leading singular behaviour of the area moments which are, in turn, obtainable by a numerical 
analysis, extrapolating finite-size data. 
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3 (/-Algebraic functional equations 



Due to the difficulty of the SAP problem, simple subclasses of SAPs have been analysed in the 
past; for example rectangles, Ferrers diagrams, stacks, staircase polygons, bar-graph polygons, 
convex polygons and column-convex polygons, among others USB E2] • All these models share 
the special property that they are directed, i.e. they can be thought of as being built up by 
adding consecutive layers of squares. This results in a special recursion for the coefficients p m ,n 
or a specific functional equation for the generating function G(x, q). Over the years it has become 
clear that the underlying structure is a q-algebraic functional equation of the form 

P(G(x, q), G(qx, <?),..., G(q N x, q),x, q) = 0, (3.1) 

where P(yo, yi, ■ ■ ■ , Vn, x, q) is a polynomial in yo, yi, . . . , yw, x and q. All exactly solved polygon 
models satisfy such a functional equation (with x conjugate to the (horizontal) perimeter and q 
conjugate to the area). For the example of bar-graph polygons, see the next section. No such 
recursion is known for SAPs, however. The limit q — > 1 in 13.1(1 generally leads to an algebraic 
differential equation, i.e. an algebraic equation in x, G(x, 1) and its derivatives. The exactly solved 
polygon models, however, result in a purely algebraic equation, i.e. in the non-trivial polynomial 
identity 

P(G(x,l),G(x,l),...,G(x,l),x,l) = 0. (3.2) 

Not much is known about the singular behaviour of general q-algebraic functional equations. 
The q-linear case, however, has been treated in some detail, see [35] and references therein. The 
polygon models of rectangles, Ferrers diagrams, and stacks are q- linear |31) . whereas the other 
exactly solved models cited above satisfy (/-quadratic functional equations. Using (|3.2() . we see 
that the perimeter generating function G(x, 1) displays an algebraic singularity at some point x c . 
The most direct way to extract the scaling function about x = x c and q = 1, whose existence we 
assume, from the defining q-functional equation, is to insert the Ansatz [12.2(1 into the functional 
equation, introduce the scaling variable s = (x c — x)(l — q)^^ and expand the functional equation 
to leading order in e = 1 — q. This will, to lowest order in e, result in a differential equation for the 
scaling function F(s). Above method is also called method of dominant balance |30|. Note first 
that the expansion of the generating function with g-shifted argument is given to leading order by 

G(q k x, q) = G (re9) (x, q) + e e F{s) + kx c e e+{1 ~^ ' F' \s) + 0(e e+2 ^-^). (3.3) 

We will analyse the case of a square-root singularity of the perimeter generating function in detail. 
(The case of a simple pole is discussed in [37]). It imposes the restrictions 




0, /0, (d x P)^0 (3.4) 
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on the functional equation, evaluated at (x, q) — (x c , 1). We introduced the abbreviations dk = d Vk 
and dk t i = d Vk d Vl for < k,l < N . The assumption of a square-root singularity also implies <fi = 29 
for the critical exponents, as follows from the relation 9 + 07 = 0. The functional equation then 
has an expansion of the form 



^ (26 \ fe d ^ + ^ +(1 - 0) *c (£ kd^ F'(s) - (d x P) s 



(3.5) 



to lowest orders in e. In order to obtain a non-trivial equation for the scaling functions, we demand 
equality of the coefficients. This results in exponents 9 = 1/3 and <j> = 2/3, and we get the Riccati 
equation 

F(s) 2 - AhF'{s) - f 2 s = 0, (3.6) 
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where the coefficients fo and f\ are given by 



The equation, together with the prescribed asymptotic behaviour, has the unique solution 

F «=- 4 4'" Ai (Gir«)' < 3 - s > 

where Ai(x) = ^ J °° cos(i 3 /3 + tx) dt is the Airy function. The scaling function is regular at the 

/ r \ 2 / 3 

origin, and its first singularity on the negative real axis is a simple pole, located at ( J s c = 
-2.338107... . 

It is important to notice that the functional form of the result is independent of the detailed 
form of the functional equation, as long as the coefficients fo and f\ are nonzero and finite! 
Thus, the scaling behaviour is determined by the type of singularity of the perimeter generating 
function. This remarkable kind of universality allows us to conjecture the functional form of the 
scaling function for self-avoiding polygons: Together with the (numerical) observation that the 
perimeter generating function of the model of rooted SAPs G^ r '(x,q) — Xj^G(x,q) also has a 
square-root singularity, one might ask if it has the scaling function (|3.8|l . 



4 Bar-graph polygons 



The model of bar-graph polygons has been defined and solved exactly in PDJ. We review the deriva- 
tion and extract the scaling behaviour as outlined above. Let G(x,y,q) denote the anisotropic 
perimeter and area generating function for bar-graph polygons. The variable q is conjugate to the 
area, while the variables x and y are conjugate to the horizontal and vertical perimeters respec- 
tively. By partitioning the set of bar-graph polygons into subsets of inflated and concatenated 
ones, it is straightforward to show that G(x, y, q) satisfies a q-quadratic functional equation in the 
horizontal perimeter variable x and the area variable q 

G{x, y, q) = G(qx, y, q) qx G(x, y,q) +y G{qx, y, q) qx + y G(qx, y, q) + qx G(x, y,q)+y qx. (4.1) 

This is depicted in Figure [3 The g-quadratic functional equation can - similarly to the treatment 




+ 



Figure 2: Diagrammatic representation of the bar-graph polygons functional equation |80) . 
of a Riccati differential equation - be linearised by means of the transformation 

G{x, y,q) = — [ —, r- - (1 + qx) , 

qx \ H{x,y,q) J 



(4.2) 



and the resulting g-linear functional equation of order two for H(x, y, q) can be solved by iteration. 
This leads to 

^ ( _ g s(l_ y)) » g (5) 
H{x,y,q)=2_^ — /„.„s > ( 4 - 3 ) 



n=0 



where {t;q) n = Ilfc=o(l — tq k ) denotes the g-product. The appearance of a g-series is a generic 
feature of solutions of q-algebraic functional equations, making it difficult to extract the asymptotic 
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behaviour for q — > 1 directly from an exact solution. For staircase polygons, the first exactly 
solved model with a q-quadratic functional equation [21 I30| , this has been done by use of saddle 
point techniques |29|. and this method can also be applied to the above solution. A more direct 
approach consists in deriving the (anisotropic) scaling function directly from the defining functional 
equation by use of the above methods. To this end, note that the defining polynomial is of the 
form P(zo, Zi,x, q) = qxz\ZQ + (1 + qx)yzi + {qx — 1)zq + qxy. The perimeter generating function 
may be found by setting q = 1 and solving the functional equation for G(x, 1). It has a square-root 
singularity at x c (y) = ■ An analysis along the above lines yields 

ti = {1 ~ y) llw m ' 5 C + »-('-»)*■(»». <«) 

and the scaling function is the logarithmic derivative of an Airy function, given by (|3.8(l . The 
critical exponents are 9 = 1/3 and <fi = 2/3. 

5 Self-avoiding polygons 

The model of rooted SAPs, G^(x, q) = x-£;G(x, q), has been analysed numerically some time ago, 
where it was found to exhibit the same exponents 9 = 1/3 and ^ = 2/3 as bar-graph (staircase) 
polygons which, in particular, implies a square-root singularity of the perimeter generating function 
(see also ^2)- With the improved series data now at hand and the refined techniques of analysis, 
we were able to confirm our conjecture about the scaling function to high precision using data 
obtained from exact enumeration for small values of perimeter, on the square lattice and on the 
triangular lattice |36| . In addition, it allowed us to analyse the first two corrections to scaling 
|37j . Here, we present results obtained from an analysis of extended enumeration data for square, 
hexagonal and triangular lattices and in addition we present data for the square lattice obtained 
from Monte-Carlo simulations. The underlying assumption is that the scaling behaviour of rooted 
SAPs may be arbitrarily well approximated by the scaling behaviour of solutions to a g-algebraic 
functional equation with a square-root singularity. The scaling function (|3.8|l implies certain 
relations between the amplitudes fk of the leading singularity (|2.5|l of the area moment functions 
gk(x). Using the Riccati equation (|3.6|l for the scaling function, it can be shown that /„ = 
Cn/f/g - ™, where the numbers c„ are given by the quadratic recursion 

_^ n— 1 

c„ + (3n - 4)c„_i + - ^ c n~rC r = 0, c = 1. (5.1) 

r=l 

The first few values are a = 1, c 2 = -5/2, c 3 = 15, c 4 = -1105/8, c 5 = 1695, c 6 = -414125/16, 
c 7 = 472200, c 8 = -1282031525/128, c 9 = 242183775, c 10 = -1683480621875/256. The coef- 
ficients fk may be estimated from series enumeration data: Assuming an asymptotic growth of 
coefficients of the area moment functions gk (x) of rooted SAPs of the form 

[x rn }g k (x) ~ {-l) k E k x- m m~ tk - 1 (m -» oo) (5.2) 

implies a singular behaviour of gk(x) of the form (|2.5|l with fk = (— l) k ^xJ k T(^k) and 7^ = 
(3k — l)/2. The constant a is defined such that p m ,n is nonzero only if m is divisible by a. Thus 
a = 2 for the square and hexagonal lattices and a = 1 for the triangular lattice. Note that the 
amplitudes E kl defined by (|5.2|) . take the same values for rooted SAPs and for (ordinary) SAPs. 
Extrapolation techniques allow us to estimate x c , 7^ and E k to high precision. This route was 
followed in [j23 EZ] where we used series data from exact enumeration up to perimeter 86 on the 
square lattice and up to perimeter 29 on the triangular lattice, where the same behaviour ()5.2JI is 
expected due to universality. 

Since then the enumeration of area-moments for square lattice SAPs has been extended to 
perimeter 100 (110) JSJj an d series have been derived for the hexagonal lattice to perimeter 140 
(156) and the triangular lattice to perimeter 58 (60) ^B], where the numbers in parenthesis refer 
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Amplitude 


Exact value 






Square 








Hexagonal 


Triangular 


Eq 


unknown 








0.56230130(2) 






1.27192995(10) 


0.2639393(2) 


E 1 


0.7957747 


X 


10" 


-1 


0.795773(2) 


X 


10" 


-1 


0.795779(5) x 10" 1 


0.795765(10) x 10" 1 


E2E0 


0.3359535 


X 


10" 


-2 


0.335952(2) 


X 


10" 


-2 


0.335957(6) x 10~ 2 


0.335947(5) x 10~ 2 


EzEl 


0.1002537 


X 


10" 


-3 


0.100253(1) 


X 


10" 


-3 


0.100255(3) x 10" 3 


0.100251(4) x 10" 3 


-B4-E0 


0.2375534 


X 


10" 


-5 


0.237552(2) 


X 


10" 


-5 


0.237557(7) x 10" 5 


0.237547(6) x 10~ 5 


E$Eq 


0.4757383 


X 


10" 


-7 


0.475736(3) 


X 


10" 


-7 


0.475749(10) x 10~ 7 


0.475724(15) x 10" 7 


E&Eq 


0.8366302 


X 


10" 


-9 


0.836624(5) 


X 


10" 


-9 


0.836652(10) x 10~ 9 


0.83660(2) x 10" 9 


E~;E% 


0.1325148 


X 


10" 


-10 


0.132514(2) 


X 


10" 


-10 


0.132519(5) x 10" 10 


0.132511(5) x 10~ 10 


EbE 


0.1924196 


X 


10" 


-12 


0.192418(2) 


X 


10" 


-12 


0.192426(8) x 10~ 12 


0.192419(8) x 10~ 12 


EgEu 


0.2594656 


X 


10" 


-14 


0.259464(2) 


X 


10" 


-14 


0.259472(12) x 10~ 14 


0.25948(4) x 10~ 14 


EiqEq 


0.3280633 


X 


10" 


-Hi 


0.328062(4) 


X 


10" 


-Hi 


0.328051(15) x 10" 16 


0.32812(5) x 10" 16 



Table 1: Predicted exact values for universal amplitude combinations and estimates from enu- 
meration data for square, hexagonal and triangular lattice polygons. 



to slightly longer scries calculated for the perimeter generating function. Details of analysis can 
be found in |171 114) and a brief summary will suffice here. Estimates for the critical points and 
exponents are obtained using the numerical technique of differential approximants |12| . For the 
square and hexagonal lattices, where a = 2 and only coefficients with m even are non-zero, we 
actually analyse the functions h k (y) = J2 m >o n nk P2m+j,ny m , where j = 4 and 6 is the perimeter 
of the smallest SAP on the square and hexagonal lattices, respectively. A singularity at x = x c 
in g k (x) thus becomes a singularity at y = x\ in hk{y)- The generating functions for the square 
and triangular lattices show no convincing signs of singularities other than the one at x c . In these 
cases we fit the coefficients to the assumed form 

Y / nk Pm,n~x- m m^~ 1 kl[E k +Y,a t /m 1+l / 2 } (m — 00). (5.3) 

n i>0 

In the hexagonal case hk{y) has an additional singularity on the negative axis at y = — x 2 _ = 
—0.412305(5) with exponents which appear to equal the exponents at y — x\ = 0.2928932..., 
as first noted in [10) . Since X- > x c the singularity at y = —x 2 _ is exponentially suppressed as 
771 — ► 00 , however for small values of m it still has a significant influence on the coefficients in the 
generating functions and we have to add an extra sequence to the fit similar to 1)5. 3[l : 



{-l) m l 2 \x-\- m m^- 1 k\[E^+^b l /77i 1 



i>0 



We obtain several data sets by varying the number of terms used in the fit. Only the even 
coefficients are used in these fits for square and hexagonal SAPs. To obtain the final estimates we 
do a simple linear regression on the data for the amplitudes as a function of 1 jm extrapolating to 
1/m — > 0. We estimate the error in the amplitude from the spread among the different data sets. 
In this way, we obtain the results for the amplitude combinations listed in Table ^ which have 
been shown to be independent of the underlying lattice [HI 136) . The scaling function prediction 
for these numbers is 



p p 2fc-l _ c 2fc (3fc-2)! 2 fc _ C 2fc+1 - 

EikE Q ~- 4n3k (6A ._ 3)! . E 2k+1 E ~ ( 3 fe)!7r 3 fc+ i 2 6fc+2 ( i ' 4 ' 



for fceN, where we used the known result that Ei = J- 0], and the numbers c k are defined in 
l)5.1|) . It is clear that the estimates for the first 10 area weighted moments are in perfect agreement 
with the predicted exact values. 

We also extended our previous approach by analysing SAPs up to perimeter 2048 by a Monte- 
Carlo simulation. An effective method of Monte-Carlo SAP generation of a given perimeter m has 
been described in ■ Each MC step consists of an inversion or a certain reflection of a randomly 
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chosen part of the SAP. Polygons which are no longer self-avoiding are rejected. The check for 
self- avoidance (by marking the corresponding lattice points) as well as the area determination 
("on the fly") can be done in time proportional to the polygon length. For a given perimeter m, 
we took a sample of 10 6 polygons, with at least 10 x m MC update moves between consecutive 
measurements. This leads to an estimated statistical error of 1 part in 1000. For perimeter 2048, 
this took two days of CPU time on a 1.5 GHz PC. We used the random number generator ran2 
described in |34| . From the Monte-Carlo measurements of polygons of a given perimeter m, we 
extracted ratios Dk / D k of the mean area moment amplitudes Dk defined by 



fe\ _ 2^2n nk Pm.n 



Dkklm 



2l,h 



(m — » oo), 



(5.5) 



where p m ,n is the number of sampled polygons of area n, and v = 3/4 follows from (|5.2|l . Note 
that by using MC sampling we cannot measure the coefficients Ek of (|5.2() directly but only the 
ratios D k = E k /E . We tested the accuracy of the random number generator by comparing values 
for the area moments obtained by simulation with their exact values, which are available up to 
perimeter 100 from exact enumeration. Assuming the above scaling function results, for k G N, in 
the values for the amplitude ratios 



D 2 k/D\ k = -c 2k 



(3fc-2)!2 4fe - 2 



(6k-3)lir k 



D 



2fc+l 



/ n 2fc+l _ 



1 



(3fc)!2 2fc 7r fe ' 



(5.6) 



where the numbers c k are defined in i|5.1|) . We computed the values (a ) m /(fc!(a)JjJ for m = 
64, 96, 128, 192, 256, 512, 1024, 2048, and a typical data plot is shown in Figure El 




0.527 



0.526 



0.0025 0.005 0.0075 0.01 0.0125 0.015 



Figure 3: MC estimates of (a 2 ) m /(2(a)^) plotted against 1/m. The estimated statistical relative error 
is one part in 1000. 



We then extrapolated the numbers Dk/D k by a least squares fit to (a ) m /(kl{a) k n ) ~ Dk/D k + 
am -1 +C2m -3 / 2 asm-t oo, in analogy with the series extrapolation for exact enumeration data. 
The corresponding results for the first ten amplitude combinations are given in Table together 
with their theoretical values. The relative error is less than 1 part in 1000. 

The value of D\ is estimated to be D\ as 0.141776. We can compare this to our previous analysis 
by noting that D\ = Ei/Eq w 0.14152105(1), where Ex = j- has been derived by conformal field 
theory arguments 0, and E = 0.56230130(2) together with x c = 0.379052277757(5) has been 

determined in [171111 E3 Noting that f = -2^y^ and f x = -x c ^- = -^we arrive at the 

conjectured form of the scaling function for rooted SAPs 



F«(s) 



ds 



logAi — {2E y s 



(5.7) 
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Amplitude 


Exact value 




MC data 


D 2 /D'{ 


0.530516 Xl0" 


-0 


0.530526 xlO" 


D 3 /Dl 
Di/Df 
D 5 /Dl 


0.198944x10" 


-0 


0.198954x10- 


0.592380x10" 


-1 


0.592445x10" 


0.149079x10- 


-1 


0.149105x10- 


De/Df 


0.329452x10- 


-2 


0.329534x10- 


Eh/Dl 
D s /Dl 


0.655743x10- 


-3 


0.655959x10- 


0.119654x10- 


-3 


0.119705x10- 


Dc,/D\ 


0.202754x10- 


-4 


0.202863x10- 


D w /D\° 


0.322149x10" 


-5 


0.322376x10- 



Table 2: Comparison of prediction for amplitude ratios against MC data. The estimated statistical 
relative error is one part in 1000. 

with exponents 6 = 1/3 and <fi — 2/3. (Formula (22) in |36| is correct up to a minus sign.) The 
conjectured form of the scaling function is then obtained by integration and is 




with exponents 8 — 1 and <f> = 2/3. The parameters for the hexagonal lattice are a = 2 and 
x c = 1/ *J~(2~-\- \/2) (known exactly from the work of Nienhuis 26 ) and for the triangular lattice 
a = 1 and x c = 0.2409175745(3). 

6 Conclusion 

We discussed the scaling behaviour of g-algebraic functional equations which underlie the critical 
behaviour of exactly solvable planar polygon models, counted by perimeter and area. This led to 
a prediction of the scaling function of the unsolved model of self-avoiding polygons on the square 
lattice, which we verified to numerical precision. The prediction also implies certain amplitude 
combinations for all area moments. It was previously argued [5] that these combinations should 
be universal, i.e., independent of the underlying lattice. Our investigations on the triangular and 
on the hexagonal lattices support the universality hypothesis, see also It is possible to extend 
the above analysis to include the first few corrections to scaling |37) . 

The question arises whether the SAP area statistics can be found in other models as well. One 
candidate is the hull of planar Brownian motion, which was observed in 1984 to have the same 
fractal dimension as SAPs |25| . Interestingly, this observation has been proved recently using 
methods of stochastic processes [20] ■ It might thus be possible to derive the area distribution 
rigorously for this model, as well as a corresponding prediction for the SAP area distribution. 

Since g-algebraic functional equations are the underlying mathematical structure of (exactly 
solvable) planar polygon models counted by perimeter and area, they may serve to classify possible 
scaling behaviour of these polygon models, g-algebraic functional equations exhibit various types 
of scaling behaviour, depending on the specific type of singularity of the perimeter and area 
generating function |37j . Since the techniques applied so far only make use of formal power series 
expansions, one might ask for a proof of existence of a uniform asymptotic expansion. Details of a 
proof may rely on techniques used in the related problem of the behaviour of differential equations 
containing a parameter [22] about singular points. Recently, the scaling function prediction has 
been rederived using field-theoretic methods and generalised to higher order critical points 0. It 
would be interesting to see whether the proposed scaling functions of arise in this framework 
as well. It may be possible that g-algebraic functional equations also provide insight into models 
of polygons with interaction, q then being the activity of interaction. There is at least one such 
exactly solvable model: self-interacting partially directed walks [2"%| . 
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Another possible application might be to models of cluster hulls, which appear in percolation 
or in spin models as boundaries of spin clusters in the Ising model or in the Potts model, as 
described in [7]- 
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